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Direct observation of the tube model in F-actin solutions: 
Tube dimensions and curvatures 
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Abstract. - Mutual uncrossability of polymers generates topological constraints on their confor- 
mations and dynamics, which are generally described using the tube model. We imaged confine- 
ment tubes for individual polymers within a F-actin solution by sampling over many successive 
micrographs of fluorescently labeled probe filaments. The resulting average tube width shows the 
predicted scaling behavior. Unexpectedly, we found an exponential distribution of tube curvatures 
which is attributed to transient entropic trapping in network void spaces. 
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Introduction. — One of the central questions of poly- 
iper physics is how the properties of polymeric materials, 
such as fluids or gels, derive from the characteristics of 
lihe individual macromolecules forming the material. In 
polymer solutions, entanglements between neighboring fil- 
aments impose topological constraints on the Brownian 
motion of individual polymers. These constraints, how- 
ever, are notoriously difficult to incorporate into a theo- 
retical description of the system. Arguably the most suc- 
cessful phenomenological approach to this problem is the 
tube model introduced by Edwards [1] . It reduces the com- 
plex many-filamcnt problem to a tractable single-filament 
problem, where the topological constraints imposed by the 
surrounding polymers are condensed to the concept of an 
impenetrable tube. This assumption is at the heart of 
most theories of macromolecular behavior such as the rep- 
tation model [2] . For solutions of linear polymers that are 
not cross-linked the confinement tube is a dynamic entity 
that is defined on intermediate time scales ranging from 
the moment of first contact between neighboring polymers 
to the characteristic time for tube remodelling. However, 
this time window can be of substantial size [2]. 

Up to now the tube model has remained mostly a con- 
ceptual device for calculations with few direct experimen- 



tal observations [3,4]. Due to this sparseness of data 
most theories on collective properties of semi-dilute or 
dense polymer solutions rely on the assumption that fila- 
ments are enclosed in confinement tubes with fixed widths 
throughout the sample. Moreover, the ensemble of con- 
finement tubes is usually assumed to qualitatively show 
the same conformation statistics as an individual polymer 
in free solution. In the present study, we challenged these 
assumptions by directly imaging confinement tubes of flu- 
orescently labeled actin filaments embedded in semi-dilute 
solutions of unlabeled filaments. 

Actin filaments (F-actin) form by self-aggregation of the 
protein actin and exhibit lengths of some tens of microns. 
In solution, these filaments undergo thermal shape fluc- 
tuations but, due to their bending stiffness, maintain an 
average direction and are therefore referred to as semi- 
flexible; their persistence length ranges from 15-20 fim 
for phalloidin stabilized filaments [5-7] to approximately 
10 /im for unstabilized filaments [8]. At present, the prop- 
erties of isolated filaments are reasonably well understood 
in terms of the wormlikc chain model [9] but the relation 
between microstructure and viscoelastic properties of en- 
tangled solutions is still a challenging problem [10, 11]. 

For semiflexible polymers, a typical tube is pictured as a 
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thin and relatively straight cylindrical pore. A thermally 
fluctuating filament confined within such a tube will en- 
counter collisions with the tube wall, whose average dis- 
tance is given by Odijk's deflection length [12]. Expanding 
this concept Scmcnov derived a scaling law for the depen- 
dence of tube width, w, on monomer concentration, c [13]. 
Basically, he exploited the fact that fluctuating filaments 
explore non-overlapping regions of space and predicted 

w oc c"3/^ . (1) 

In contrast, the mean distance between neighboring fila- 
ments, .J, i.e. the mesh size, scales for purely geometrical 
reasons like c~^/^. Thus mesh size and confinement tube 
width are genuinely different length scales. Confinement is 
usually described by a harmonic tube potential which has 
been confirmed by experiments [14] and computer simula- 
tions [15]. 

In this Letter, we present a quantitative study of con- 
finement tubes in entangled F-actin networks. We fluo- 
rescently labeled a trace amount of actin filaments with 
rhodamine-phalloidin, embedded those in a solution of 
unlabeled filaments and observed individual filaments by 
confocal laser scanning microscopy. A typical micrograph 
is shown in Fig. la. Upon averaging time series of such 
images, we obtained an intensity distribution representing 
the probability density of a filament within its tube (Fig. 
lb). This approach allowed us to directly characterize the 
size and conformation of individual confinement tubes in 
real space at optical resolution. 

Experimental. — Actin was prepared as described in 
the literature [16]. Monomeric actin was kept in G-buffer 
(2.0 mM Tris-HCl, pH 8.0, 0.2 mM CaCla, 0.2 niM ATP, 
0.2 mM DTT) at 4°C and used for experiments within 
72 h. Polymerization was initiated by adding 1/10 of the 
sample volume of 10 x F-buffer (20 mM Tris-HCl, pH 7.5, 
2 mM CaClz, 5 mM ATP, 2 mM DTT, 20 mM MgClg, 1 M 
KCl). The samples were then incubated for 2 h at 37° C. A 
part of the F-actin was labeled with rhodamine-phalloidin 
(TRITC-phalloidin) at a molar ratio of 1:1. Labeled and 
unlabeled F-actin solutions of the same concentration were 
then gently mixed at 1:1000 ratio for 15 min in a slow 
(few rpm) rotary shaker. The resulting solution was gently 
pipetted onto a micro slide and overlayed with a cover 
slide which was held at a distance of 100 /zm by a spacer 
of adhesive tape. Experiments were performed after an 
equilibration period of 15 min. 

Filaments were observed using an inverse confocal mi- 
croscope (LSM510 equipped with a water immersion ob- 
jective, C-Apochromat 40/1. 2W corr, Carl Zeiss, Jena, 
Germany) . Excitation wavelength was 543 nm. Only fila- 
ments residing in a plane perpendicular to the optical axis 
at a distance of 40-60 /im from the cover slide were con- 
sidered. The thickness of the optical slice was typically 
set to 2.8 /im. A field of view containing one labeled fil- 
ament was scanned within 0.1-0.3 s. Images were taken 
at 1 s intervals within a time period of 2-3 minutes. On 
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Fig. 1: (a) Micrograph of a fluorescently labeled F-actin within 
a solution of concentration c =14.5 /iM. (b) Average of 100 
successive images taken within 126 s. Overlayed is the center 
line of the averaged contour (black line). Scale bar 5 /^m. (c) 
The width of the averaged intensity distribution corrected for 
optical resolution effects. Arc length runs from left to right in 
the micrographs. 

this time scale we observed only transverse filament fluc- 
tuations. Time-lapse sequences of about 25 filaments were 
collected for each concentration. 

As contour characterizing a confinement tube we used 
the center line of the intensity distribution in the time 
averaged micrographs, cf. Fig. lb. We determined this 
contour by the following algorithm: The image was bina- 
rized with a global threshold. The resulting areas were 
skeletonized to obtain a first estimate of the contour at 
pixel accuracy. Then the end points of this line were kept 
fixed and an active contour algorithm based on energy- 
minimization was applied on the time averaged micro- 
graphs as described in Ref. [17]. This resulted in a smooth 
line representing the contour of the confinement tube; see 
Fig. lb. All image processing routines were developed us- 
ing MatLab (Release 14, The Mathworks, Natick, MA). 
The width of the confinement tube, (Joxp, was determined 
as the standard deviation of a Gaussian fit to the intensity 
distribution of the averaged image along lines perpendic- 
ular to the contour. A result from this procedure is shown 
in Fig. Ic. 
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Fig. 2; Dependence of the average tube widths, w, on actin 
concentration, c. Symbols: measured points, error bars: stan- 
dard deviations of the mean tube widths of the filaments, full 
line: best fit to a power law (resulting exponent -0.61), dashed 
line: best fit with exponent -1/2. 



Experimental Results. — Our first objective was 
to analyze the functional dependence of the average tube 
width w on the actin monomer concentration c. The 
data for the width CTexp were corrected for the optical 
resolution^ clsf to give an undistorted tube width w = 
('''cxp ~ '''lsf)^^^- Wc also took into account that a so- 
called critical concentration (1.9 /iM) of actin monomers 
forms a pool of unpolymcrizcd molecules in equilibrium 
with the actin filaments, i.e. in Eq. 1 the concentration of 
actin within filaments was taken to be 1.9 fiM less than 
the total actin concentration. Average widths of the con- 
finement tubes were 0.34, 0.26, 0.22, and 0.18 /im at total 
actin concentrations of 14.5, 21.8, 29.1, and 36.3 fiM, re- 
spectively. These data conform quite well to a power law 
curve with a best fit for an exponent of -0.61, cf. Fig. 2. 
An exponent -1/2, as expected for the mesh size, is clearly 
not compatible with our data. To our knowledge, this 
result is the first experimental confirmation of the pre- 
dicted Semenov scaling exponent of -3/5, Eq. 1, showing 
that indeed tube width and mesh size are distinct length 
scales [13]. Moreover, even the absolute values of the tube 
width as predicted by recent theoretical studies [15, 18] 
are of comparable but lower magnitude. This is consis- 
tent since these treatments do not include collective fluc- 
tuations of the entangled solution and therefore give lower 
bounds. 

So far our results nicely confirmed previous theoretical 
predictions on the average tube width. However, contrary 
to common belief, fluctuations of the tube widths were 
pronounced with standard deviations of Ri35%; cf. Fig. Ic. 
The length scale over which tube width fluctuations oc- 
curred was studied via the autocorrelation function of 



^The optical resolution was calculated by convoluting the point 
spread function measured on 100 nm fluorescent beads with a line. 
The resulting line spread function (LSF) was well approximated by 
a Gaussian of width ctlsf, which on average was 0.18 fira 
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Fig. 3: Distribution functions of the curvatures of the confine- 
ment tubes. -I- 36.3 /iM, x 29.1 fxM, o 21.8 fiM, and o 14.5 nM. 
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defined as ACFU) - 
this equation the subscript "exp" was dropped for the sake 
of brevity. The initial decay parameters of the autocorre- 
lation functions of Ccxp were 0.33, 0.49, 0.60, 0.68 ^m~^ 
at concentrations of 14.5, 21.8, 29.1, and 36.3 /xM, respec- 
tively. Thus tube width fluctuations occurred on length 
scales of few microns. The faster decay of correlations in 
the tube width at higher concentrations is a consequence 
of the higher obstacle density. 

To go beyond a characterization of a tube in terms of its 
typical size we investigated tube conformations. Specifi- 
cally, we asked for the probability distribution function 
of tube curvatures. In the absence of excluded volume it 
seems obvious that the probability distribution of poly- 
mer conformations in solution is identical to that of free 
filaments. Sometimes, it is argued that even the (coarse- 
grained) conformations of the tube contour should follow 
such a free distribution [19]. Since the effective free en- 
ergy of the wormlike chain model is quadratic in the local 
curvature this would imply a Gaussian distribution for the 
curvatures of confinement tubes. To our great surprise we 
found a curvature distribution with a pronounced expo- 
nential tail; see Figs. 3 and 4 top. 

The initial decay length of the curvature distribution 
shows little to no dependence on concentration; the decay 
lengths of the distributions were 8.0, 8.0, 8.2, and 9.9 fim 
at the concentrations of 36.3, 29.1, 21.8, and 14.5 /iM, 
respectively. This is to be expected since small curvatures 
are dominated by bending stiffness. 

An exponential distribution decays much slower towards 
large values than a Gaussian. This is reflected in our fre- 
quent observation of highly bent fllaments. For isolated 
filaments, similarly bent contours are extremely rare and 
were never observed [5-7]. What is causing such a quali- 
tative change in chain conformations? At this point sev- 
eral reasons for the occurrence of highly bent shapes, be- 
sides topological interactions, were conceivable. We could 
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rule out that filaments assumed such conformations dur- 
ing sample preparation and remained kinetically trapped 
by additional experiments where the time delay between 
successive micrographs of actin filaments was set in the 
range from 10 to 30 sec. With these longer time delays 
we were able to track individual filaments for up to one 
hour. On these time scales we observed several initially 
straight filaments going into highly bent shapes. This ob- 
servation indicates that the highly bent shapes causing the 
exponential decay of the distribution of curvatures of the 
confinement tubes are an intrinsic property of the system 
and not just kinetically trapped remainders of the prepa- 
ration procedure. 

Simulation. — We resorted to Monte Carlo simula- 
tions to elucidate the underlying causes for the observed 
exponential distribution of tube curvatures. We con- 
sidered an initially equilibrated test polymer in a two- 
dimensional observation plane surrounded by randomly 
placed point-like obstacles representing the intersection 
points of neighboring polymers with that very plane. The 
tube contour was obtained by observing the test polymer's 
fluctuations in the presence of harmonically fluctuating 
point obstacles. Since these obstacles represent themselves 
polymers with identical properties, we tuned the fluctua- 
tion amplitudes of the obstacles to equal those of the test 
polymer to obtain a fully self-consistent model. Sampling 
over filaments in different network configurations and av- 
eraging over all statistically allowed conformations in a 
given tube resulted in a curvature distribution that also 
featured a pronounced exponential tail. As this tail shows 
no significant dependence on concentration in the experi- 
mental parameter range we compare the simulation results 
to the combined experimental data for all concentrations 
and find good agreement; see Fig. 4 top. 

We therefore concluded that the observed exponential 
distribution of tube curvatures is indeed characteristic for 
semi-dilute solutions of semifiexible polymers. Compared 
to free filaments, both smaller and larger curvatures are 
more frequent. An increased frequency of small curva- 
tures is due to the coarse-graining procedure involved in 
determining the tube contour: averaging over all topo- 
logically allowed polymer conformations within the tube 
(simulations) or performing a time average (experiment) 
results in a smoother contour and shifts weight of the 
probability distribution towards small curvatures. The 
origin for the relatively large number of highly bent fil- 
aments and concomitantly the exponential tail of the cur- 
vature distribution is far less obvious. To avoid compli- 
cations from coarse-graining we analyzed the distribution 
function of the confined filament itself. Even though not 
easily accessible to experiments like the tube contour, it 
allows a direct comparison to free filaments. Again, a 
fully self-consistent simulation of the polymer configura- 
tion resulted in an exponential tail in the curvature dis- 
tribution function, cf. Fig. 4 bottom, while small curva- 
tures remained unaffected. Monitoring polymer segments 
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Fig. 4: Top: Curvature distribution of confinement tubes of 
probe filaments in a dynamic network (experimental data and 
simulations) and distribution function of a free polymer (simu- 
lations). Bottom: Simulated curvature distributions of free fila- 
ments, filaments within a self-consistently fluctuating network, 
within a network of flxed and a network of strongly fluctuating 
obstacles. 

with unusual high curvatures during the simulation en- 
abled us to attribute this effect to network environments 
where a favorably bent test polymer can protrude into a 
large void space thereby realizing a higher entropy^. Upon 
immobilizing the obstacles the effect is enhanced, while for 
strongly fiuctuating obstacles it decreases, cf. Fig. 4 bot- 
tom. This is due to a smaller effective void volume and 
is expected to recover the free polymer case when strong 
obstacle delocalization renders the network comparable to 
a gas. The phenomenon bears some similarity with "en- 
tropic trapping" [20, 21] observed for flexible polymers in 
random environments. 

Discussion. — Clearly, in the absence of excluded vol- 
ume, time averaging of a single test polymer over infinite 
time as well as averaging over an ensemble of polymers at 
a given time both have to reproduce the Gaussian distri- 
bution of free filaments. The fact that our results do not 
reflect this thermodynamic equilibrium might be disturb- 
ing at first glance, but on closer inspection turns out to 
be a consequence of the tube concept and the time scales 
involved. 

It is important to realize that for a solution of not 
cross-linked polymers the tube model itself is a concept 
to describe polymer dynamics and not equilibrium poly- 
mer conformations. In polymer solutions the confinement 
tube of a given polymer is defined as the space explored 
by the polymer's fiuctuations on intermediate time scales 
ranging from the moment when the polymer first experi- 
ences topological restrictions by nearest neighbors to the 
time scale of tube remodelling. Both points in time can be 
estimated from experiment. In dynamic light scattering a 

^H. Hinsch and E. Frey, to be published 
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cross-over from single filament to restricted dynamics is 
seen at about 10 msec [22]. The time scale of tube re- 
modelling can be estimated from experiments on actin 
rcptation by the Sackmann group [23, 24]. From their 
data we deduced rcptation times of 15 min for a 1 ^m 
filament, 9 days for a 10 jim filament and 3 years for a 
50 fim filament. As reliable curvature determination re- 
quires a polymer length of at least 8 /im, experimental 
observation of terminal relaxation by rcptation is entirely 
impossible. In essence the same holds for a molecular 
dynamics simulation (even with coarse grained models) 
of the process where simulation times by far exceed the 
possibilities of even the most recent supercomputers due 
to the enormous number of different topologies. Interest- 
ingly, for actin filaments of some micrometer length the 
often neglected treadmilling motion occurring at a rate 
of approximately 2 fim per hour [25, 26] should dominate 
terminal relaxation. Yet, even a relaxation time scale of 5 
hours for a 10 fini filament is far beyond the time window 
of light microscopic observations on such delicate samples. 
From these literature data we estimate the time window 
during which a confinement tube is well defined to range 
from 10 msec to 5 hours for not stabilized actin filaments 
of 10 fim length. The latter value increases to 9 days for 
stabilized actin networks where the treadmilling motion 
of actin is abolished by phalloidin. These values demon- 
strate that actin filaments in solution are trapped in their 
localization tubes during almost the entire experimentally 
accessible time window. Consequently, all observations on 
this time scale are crucially influenced by the tube's prop- 
erties. 

This also holds for the average performed in both exper- 
iments and simulations where we start from an ensemble 
of initial configurations of test polymers embedded in the 
network, and sample their configurations over a given fi- 
nite time window. Hereby we avoid the infeasible sampling 
of the complete phase space with one slowly reptating test 
polymer, but for each initial condition we sample only over 
the part of phase space that is represented by the associ- 
ated localization tube. The ensemble has to be sufficiently 
large to incorporate all relevant areas of the phase space. 
Each member of this ensemble may now be characterized 
in terms of its "topology" given by the topological con- 
straints imposed by its neighbors. These topologies rep- 
resent a partitioning of phase space. If the dynamics over 
the finite time window would conserve the topology, sam- 
pling again would necessarily have to give a Gaussian dis- 
tribution. Actually, however, the polymer dynamics inside 
their localization tubes as well as the Markovian dynamics 
employed in the simulations are metrically transitive, i.e. 
the partitioning is not maintained under its dynamics. In 
particular, local network configurations where large void 
spaces are explored as explained above by means of rcpta- 
tion or "breathing" result in axial motion of the polymer's 
ends and thus allow for a modification of the topological 
partitioning. Hence, entropic trapping and the concomi- 
tantly highly bent configurations of filaments and tubes 



are intimately connected to transient non-equilibrium dis- 
tribution functions that are observed on time scales well 
below large scale rcptation sets in. 

In summary, we studied the conformations of confine- 
ment tubes in entangled solutions of scmiflcxiblc polymers 
by direct visualization. Our results show that confinement 
tubes are very real objects for these systems. The depen- 
dence of the average tube width on actin concentration 
conforms to Semenov's scaling prediction [13]. In exper- 
iment and simulation we identified coinciding exponen- 
tial tails in the curvature distribution of the confinement 
tubes. This deviation is due to transient entropic trapping 
and is observed on time scales below large scale rcptation. 
These findings quantify and extend the tube concept and 
demonstrate that solutions of entangled semiflexible poly- 
mers arc not accurately described by standard concepts of 
equilibrium thermodynamics on experimentally relevant 
time scales. 
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